Anisotropic Caging of Interstitial Vortices in Superconductors with a Square Array of 

Rectangular Antidots 
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We investigate anisotropy in the vortex pinning in thin superconducting films with a square array 
of rectangular submicron holes ( "antidots"). The size of the antidots is chosen in such a way that 
it corresponds to a saturation number n s = f , i.e. each antidot can trap at most one flux quantum. 
Therefore, interstitial vortices, appearing when the magnetic field exceeds the first matching field, 
are "caged" at the interstitial positions by the repulsion from the saturated antidots. We observe 
an overall higher critical current I C (H) when it is measured parallel to the long side of the antidots 
than the I C (H) along the short side of the antidots. Although the pinning force, exerted by the 
empty antidot on the vortex, turns out to be isotropic, our I C (H) data indicate that the caging force, 
experienced by the interstitial vortices and provided by the array of saturated antidots, is strongly 
anisotropic. 



I. INTRODUCTION 

Type-II superconductors (SC's) with nano-engineered 
pinning arrays are good model systems to study the fun- 
damentals of vortex pinning, since, within the limits of 
the lithographic process used for their fabrication, they 
offer the ability to tailor the pinning potential at will. 

In SC's with a periodic pinning array (PPA), the vor- 
tices form regular geometrical patterns, commensurate 
with the pinning array, at the integer H n and fractional 
Hp/ q matching fields^. This strongly reduces the vortex 
mobility and consequently increases the critical current 
I c . These commensurability effects have been intensively 
studied for square or triangular arrays of antidotsiSi^ 
or magnetic dots£*£*£. The four-fold symmetry in the pin- 
ning properties, which is induced by a square PPA, can 
be broken by using rectangular antidots. When no inter- 
stitial vortices are present, the critical current I C {H) is 
higher in these systems when the current is applied paral- 
lel to the long side of the antidots, than when it is applied 
along their short side, due to a pronounced anisotropy in 
the vortex-vortex interaction^. 

The size of the antidots determines the saturation 
number n s that indicates the maximum flux n<fio that 
can still be trapped by the antidot. When the number 
n of flux quanta applied per antidot crosses over from 
n < n s to n > n s , the pinning potential of the antidots 
changes from attractive to repulsive. 

In the present work, we have made a film with an ar- 
ray of antidots with n s = 1. This implies that interstitial 
vortices, caged by the resultant repulsion of the saturated 
antidots, appear when the magnetic field exceeds the first 
matching field. This situation is very different from pre- 
vious workS on rectangular antidots, where no intersti- 
tial vortices were present. We investigate the anisotropy 
of the vortex-vortex interaction between interstitial vor- 
tices, and of the caging potential, felt by the interstitial 
vortices and provided by the repulsive potential of the 



saturated antidots. 



II. EXPERIMENTAL DETAILS 

We patterned a SC Pb film in a 5 x 5 mm 2 cross- 
shaped geometry (see Fig. to allow electrical trans- 
port measurements in two perpendicular current direc- 
tions. The central part of the cross consists of two 
300 /im wide strips containing the square array of rect- 
angular antidots (see dark gray area in Fig. " a)). In 
both strips, the long side of the antidots points in the 
y-direction. This pattern was prepared by electron- 
beam lithography in a polymethyl metacrylate/mcthyl 
metacrylate (PMMA/MMA) resist bilayer covering the 
Si0 2 substrate. A Ge(10 A)/Pb(500 A)/Ge(200 A) film 
was then electron-beam evaporated onto this mask while 
keeping the substrate at 77 K. A liftoff process in warm 
acetone finalized the preparation. Figure ^c) shows an 
atomic force microscopy (AFM) topograph of a 6 x 6 /im 2 
area of the square antidot lattice with a period of 1.5 /im. 
The antidots have a rectangular shape (0.5 x 0.65 fim 2 ) 
with rounded corners. As shown in the schematic rep- 
resentation of a unit cell of the array in Fig. [Jb), the 
superconducting paths between the antidots are 0.85 fjjn 
and 1 /xm wide, for the x- and y-direction, respectively. 

The transport measurements were performed in a 4 He 
cryostat with the magnetic field applied perpendicular to 
the film surface. The superconducting critical tempera- 
ture was found to be T c — 7.225 K. To determine the co- 
herence length £(0), we measured the linear T C (H) phase 
boundary of a reference film, evaporated under identical 
conditions, without any in-plane nanostructuring. From 
the T C {H) slope, anc£ 
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we find £(0) = 37 nm. Using the dirty limit (£ < £q) 
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FIG. 1: Layout of the Pb film with a square array of rectan- 
gular antidots. (a) Cross-shaped geometry of the sample to 
allow for transport measurements in the x- and y-direction. 
(b) Schematic representation of a unit cell (1.5 x 1.5 /jm 2 ) of 
the array, (c) Atomic force micrograph of a 6 x 6 fim 2 area 
of the antidot lattice. 



expression £(0) = 0.865y / £o7 and the BCS coherence 
length for Pb, £o= 83 nmi£, we determined the elastic 
mean free path I — 22 nm. We then derive the pene- 
tration depth A(0) = 47 nm from the dirty limit expres- 
sion A(0) = 0.66Al^/£o/^i using Al=37 nm as the Lon- 
don penetration depth— The presence of antidots in a 
superconducting film increases the penetration depth—, 
A(0) = , A< ~°" ) = 56 nm, where S a is the area occu- 

Y 1 — 25 a / St 

pied by the antidots and St is the total area of the film. 
The Ginzburg-Landau (GL) parameter n = A(0)/£(0) is 
therefore K rs 1.5 > l/v2 R> 0.707. We conclude that the 
film with the array of rectangular antidots is a type-II 
superconductor. 



III. RESULTS 

In Fig. we show the critical current versus field 
curves I C (H), normalized to the value at zero field, 
I co = I c (H = 0), at two temperatures (T/T c = 0.989 in 
(a) and 0.992 in (b)). We have used a voltage criterion of 
Vcrit = 100 /xV. The absolute values of the critical cur- 
rent density at zero field measured at T/T c — 0.992 is of 

the order of ~ 4 • 10 8 ^- in both directions. The field axis 

m 

is given in units of the first matching field Hi, the field 
at which the density of o - v °rtices equals the density of 
antidots: fioHi = &o/d 2 = 9.2 Oe, with d = 1.5 /im the 
period of the antidot lattice and 4>o = ^ the supercon- 
ducting flux quantum. For comparison, we have depicted 
in Fig.Efc) the I C {H) curves at T/T c = 0.992 for a super- 
conducting film with a square array (period 1.5 /im) of 
rectangular antidots with a larger size (0.6 x 1.15 /im 2 ) 8 . 
Due to the presence of the antidot latticed, both the 
lex (H) and I cy (H) data show pronounced maxima at the 
integer matching fields Hi, H 2 = 2_ff 1 , and iJ 3 = ZHi. 
Moreover, as it is most clearly seen in the curve measured 




-2-1012 



H/H, 

FIG. 2: (a) and (b) Normalized critical currents and I cy 
as a function of normalized magnetic field H/Hi for the sam- 
ple shown in Fig. Q measured with a current in the x- (open 
symbols) and the y-direction (filled symbols). Data are pre- 
sented for two temperatures: (a) T/T c = 0.992 and (b) T/T c 
= 0.989. (c) I CX (H) and I cy (H) curves at T/T c = 0.992 for a 
film with a square array of rectangular antidots with a larger 
size (0.6 x 1.15/mr 2 ). The data were obtained for H > 
(open and filled symbols) and then symmetrized for clarity 
for H < (dashed and solid lines) [adapted from Refi-]. 



with a horizontal current (open symbols in Fig.ETa-b)). 
we find a sharp drop in critical current after Hi, in- 
dicating the presence of mobile interstitial vortices for 
H > Hi 3 . This means that the antidots can trap at 
most one flux quantum, in agreement with the calcu- 
lated saturation number— n s (T) = 2 £(T) = ^-45 ~ 1 
for the temperatures under consideration, where we have 
used for a the average of the two sides of the rectangular 
antidot. 

Let us first focus on the critical current behaviour in 
the field range where no interstitial vortices are present, 
H < Hi. For these fields, the critical current I cy {H) 
(filled symbols in Fig. |2| is considerably enhanced com- 
pared to I CX (H) (open symbols) except exactly at the 
first matching field, where an identical value for I c (Hi) 
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is found for both current directions. Another striking 
feature is that rational matching peaks, although clearly 
present for I\\x (e.g. at H/Hi = 1/4 and 1/2), are al- 
most completely suppressed when 7||y. These findings 
are fully consistent with previous experiments 8 on films 
with rectangular antidots where the field range, where 
no interstitial vortices were present, extended up to Hi 
(n s -2) (see Fig. He)). 

At higher magnetic fields H > Hi, i.e. when intersti- 
tial vortices appear in the superconductor with n s = 1, 
the critical current I cy (H) (filled symbols in Fig.[2Ia-b)) 
is significantly larger than I CX (H) (open symbols) for ev- 
ery magnetic field, even at the integer matching fields 
H2 and H3. In the film with n s = 2 (Fig. Hfc)), in- 
terstitial vortices are not present at H^. Therefore, the 
critical current at the second matching field is isotropic, 
i.e. Iex{H%) = I c y(H2). Rational matching peaks are 
again only found for J||x, namely at H/H\ — 1.5 (open 
symbols in Fig. |2a-b)), whereas no rational matching is 
found for /|jy. 



IV. DISCUSSION 

The observed differences in critical current I C {H) along 
the two current directions - horizontal and vertical - are 
due to an anisotropy of the vortex mobility in the sam- 
ple along the two in-plane directions (x and y). For a 
certain Lorcntz force, this mobility depends on the pin- 
ning force experienced by the vortex, and on the vortex- 
vortex interaction force between neighbouring vortices. 
Interestingly, since the sample in the present work has 
n s = 1, the critical current probes the mobility of two 
kinds of vortices, depending on the magnetic field regime, 
i.e. of the vortices pinned at the antidots and of the 
interstitial vortices caged in between the saturated an- 
tidots. Indeed, at fields below the first matching field 
Hi, the critical current I C {H) yields information about 
the pinning force provided by the rectangular antidots 
and about the vortex- vortex interaction between vortices 
pinned at the antidots. At magnetic fields above Hi, 
the critical current is governed by the mobile intersti- 
tial vortices. Therefore, it provides information about 
the caging force produced by the saturated antidots and 
about the vortex- vortex interaction between two intersti- 
tial vortices. 

Let us discuss first the regime H < Hi, where no in- 
terstitial vortices are present: 

The value of I c (Hi) is a measure of the single site 
pinning force of the individual antidots. Indeed, at the 
first matching field, the vortex-vortex interactions in the 
highly symmetric vortex lattice (see Fig.|3Ia)) cancel out 
and depinning of the vortices occurs simultaneously once 
the Lorentz force exceeds the pinning force provided by 
the antidots. Since the value of the critical current at 
the first matching field is identical in the x- and the y- 
direction, i.e. I cx {Hi) = I cy {Hi) (Fig. we conclude 
that the pinning force exerted by the rectangular anti- 



dot on a single pinned 0o-vortex is isotropic along the 
two symmetry-axes of the rectangular antidot. Appar- 
ently, the aspect ratio of the rectangular antidots should 
be much higher in order to produce an anisotropic pin- 
ning forc o 8 i 13 . Keeping one side of the antidots fixed 
(semi-axis b = 0.25 /^m), and taking into account that 
the coherence length at the considered temperatures is 
£(T) sa 0.4 /im, one would have to construct antidots 
with an aspect ratio of ~ 8, i.e. with a long semi-axis 
a ~ 2 /an, in order to observe any anisotropy in the pin- 
ning force provided by the antidot^. The aspect ratio, 
required to obtain an anisotropic pinning force, is there- 
fore at least a factor of 6 higher than the we have used 
(see Fig. Hp))). 

The appearance of rational matching features is ev- 
idence of a strong vortex-vortex interaction in the di- 
rection of the Lorentz force. Since the I CX (H) curve 
(Lorentz force in the y-direction) shows rational match- 
ing and the I cy {H) curve does not, we conclude that the 
vortex-vortex interaction between vortices pinned at the 
antidots is stronger in the y- than in the x-direction. 

In the field regime H > Hi, interstitial vortices are 
present in our n s = 1 sample. Exactly at the second 
matching field H2 , a highly symmetric vortex configura- 
tion is formed, as is schematically shown in Fig. [3£b). 
Half of the vortices are strongly pinned at the antidots, 
the other half is caged in between the saturated antidots. 
The two vortex species form interpenetrating square lat- 
tices. When a current is applied along the x- or y- 
direction, the interstitial vortices will move first, since 
they are only weakly pinned^i^i^. Therefore, the critical 
current at H2 is a measure of the caging force provided by 
the saturated antidots and acting on the interstitial vor- 
tices. Since I cy (H2) > IcxiH^) (Fig.[2Ia-b)), the Lorentz 
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FIG. 3: Schematic presentation of the vortex configuration in 
the n s — 1 film with a square array of rectangular antidots (a) 
at Hi, (b) at Hi- Vortices are represented by current arrows. 

force needed to depin the interstitial vortices is larger in 
the x- than in the y-direction, and we conclude that the 
caging force, resulting from the repulsion of the saturated 
antidots, is anisotropic, i.e. stronger in the x- than in the 
^-direction. 

Moreover, we only find rational matching features in 
the I CX (H) curves in the field regime H > Hi. This 
indicates that the vortex-vortex interaction between two 
caged interstitial vortices is also anisotropic, i.e. stronger 
in the y-direction than in the ^-direction. 

Summarizing our observations, we have found that 
the vortex-vortex interaction between vortices pinned at 
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the antidots, as well as between interstitial vortices, is 
anisotropic, i.e. the interaction is stronger parallel to the 
long side of the rectangular antidots (y-direction) . The 
pinning force, experienced by a caged interstitial vortex 
and provided by the surrounding saturated antidots is 
anisotropic, i.e. stronger along the short side of the an- 
tidots (x-direction). Interestingly, the pinning force pro- 
vided by the rectangular antidots themselves is isotropic. 

Most of these features (except for the last, for which 
we refer to Ref£) can be understood from the geome- 
try of the antidot array. The superconducting strands 
between the rectangular antidots, where the screening 
currents of the vortices have to flow, are thinner between 
the adjacent antidots in the y-direction than in the re- 
direction (see Fig.GUa)). Therefore, we expect a stronger 
vortex-vortex interaction in the y-direction between vor- 
tices pinned at the antidots. Moreover, when interstitial 
vortices are present, it is clear from the schematic pre- 
sentation of the vortex lattice in Fig. Efb), that the cir- 
culating currents of the interstitial vortices are shielded 
more in the x-direction, which leads to a weaker vortex- 
vortex interaction between interstitial vortices in the in- 
direction. Finally, the anisotropic caging force can be 
understood if one considers that the rectangular shape of 
the antidots forces the circulating currents of the vortices 
trapped inside the antidots to adopt elongated stream- 
lines, rather than circular ones. Since the antidots are 
placed in a square array, the currents will overlap more 
in the y-direction (Fig.|3J). The caging potential barrier 
will therefore be higher in the x-direction, i.e. it will be 
more difficult for an interstitial vortex to pass in between 
saturated antidots in the x- than in the y-direction. 

In other words, due to the application of the array of 
rectangular antidots, the amount of superconducting ma- 
terial per unit length is smaller for the y-direction than 
for the x-dircction. Consequently, shielding of the mag- 
netic field is poorest in the y-direction, leading to a larger 
effective penetration depth A in this direction. Vortices 
therefore have a stronger vortex- vortex interaction in the 
y- than in the x-direction. 

It is interesting to compare these results with molec- 
ular dynamics simulations on a similar system, a rect- 
angular array of isotropic pinning sites that can trap 
at most one vorte:»i&. In this system, a different spac- 



ing in the x- and y-direction of the pinning sites induces 
both the anisotropy in vortex- vortex interaction between 
pinned vortices and in the vortex-vortex interaction be- 
tween vortices caged at interstitial positions between the 
saturated pinning sites. The calculated critical current 
versus field curves for a driving force in the two directions 
show exactly the same features as found in our I C (H) 
measurements. 



V. CONCLUSION 

We have investigated vortex pinning in a film with a 
square array of rectangular antidots, focussing on the 
behaviour of the interstitial vortices in this anisotropic 
environment. Since the antidots we have used have a sat- 
uration number n s = 1, i.e. they can trap at most one 
flux quantum, interstitial vortices appear when the mag- 
netic field exceeds the first matching field Hi. Although 
the aspect ratio of the antidots is very small, the caging 
force, provided by the saturated antidots and felt by the 
interstitial vortices, is stronger along the short side of 
the antidots. Moreover, we find a weaker vortex-vortex 
interaction between interstitial vortices in this direction. 
This leads to an overall higher critical current I C (H) with 
no rational matching features when it is measured par- 
allel to the long side of the antidots. This observation 
may be important for practical applications of films with 
a periodic pinning array, since it solves the problem of 
a reduced critical current density between the matching 
peaks in I C (H) at integer H n and fractional H p / q , thus 
leading to a substantial overall current enhancement due 
to the presence of the rectangular antidots. 
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